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1.  Introduction  and  Summary .  A  random  vector  follows  a  symmetric  multivariate 

normal  (SMN)  distribution  (Rao,  1973)  if  the  components  have  equal  means, 

equal  variances  and  equal  covariances  —  the  correlation  coefficient,  p, 

between  any  two  components  is  the  same  and  is  termed  the  intraclass,  equi-, 

uniform  (Geisser,  1963)  or  familial  (Fisher,  1925)  correlation.  Since  they 

arise  naturally  in  psychology,  education,  medicine,  genetics  etc.,  such 

models  have  received  considerable  attention,  e.g.  in  the  same  year  and  in 

the  same  journal  by  Geisser  &  Desu  (1968),  Gleser  (1968)  and  Han  (1968) 

and  very  recently  by  Mak  &  Ng  (1981).  Again  in  bio-sciences  (Press,  1981) 

many  organisms  e.g.  starfish,  octopus  etc.  exhibit  a  natural  symmetry  and 

consequently  give  rise  to  the  above  distribution.  The  importance  of  such 

models  in  many  tests  of  multivariate  analysis,  e.g.  in  MANOVA,  Profile 

Analysis,  Growth  curve  analysis  etc.  has  been  established  by  Huynh  and 

Feldt  (1970) .  Wilks  (1946)  and  Votaw  (1948)  have  considered  likelihood 

ratio  tests  for  these  models.  Tests  for  p  have  been  proposed  by  Srivastava 

(1965)  using  Roy's  union  intersection  principle,  by  Aitken,  Reinfurt  and 

Nelson  (1968)  and  by  Mak  &  Ng  using  the  likelihood  ratio  (LR)  principle, 

by  Rao  (1973)  using  the  canonical  form  of  the  SMN  distribution  and  by 

GOkhale  &  SenGupta  (1982)  using  the  approach  of  locally  most  powerful  (u^p)  tests. 

These  models  have  been  investigated  using  the  theory  of  products  of  problems 

by  Arnold  (1973)  and  in  the  general  framework  of  representation  theory  and 

invariance  by  Andersson  (1976).  Recently,  they  have  been  generalized,  e.g. 

by  Szatrowski  (1979)  to  incorporate  block  symmetries. 

A  random  vector  X  will  be  said  to  follow  a  standard  symmetric  multi¬ 
variate  normal  (SSMN)  distribution  if  it  follows  a  SMN  distribution  and 


additionally  the  components  have  zero  means  and  unit  variances.  Though 
the  literature  on  the  SMN  distribution  is  quite  extensive,  yet  little 
is  known  about  SSMN  distributions.  In  particular,  no  tests  for  p  is 
known  for  SSMN  distributions.  However,  such  distributions  can  occur  naturally 
in  various  ways  (Sampson,  1978).  One  such  situation  is  'when  there  are 
many  observations  individually  taken  at  different  times  on  each  scalar 
multivariate  normal  random  variable,  and  then  several  vector  observations 
are  taken  for  all  the  variables  together'.  Another  situation  is  'when 
there  are  many  missing  observations  for  individual  entries  in  the  vector 
sample'  and  complete  observation  vectors  are  treated  as  coming  from  the 
distribution  with  'known*  means  and  variances  found  by  using  vectors  with 
missing  entries.  Further,  in  many  practical  problems  it  is  necessary  to 
standardize  the  variables .  Even  the  sample  means  and  variances  are 
usually  employed  for  such  standardizations  and  then  the  resulting  variables 
behave  asymptotically,  by  Slutsky's  theorem,  as  standardized  random  vari¬ 
ables.  Such  standardizations  are  always  made  and  play  important  roles  in 
the  techniques  for  reduction  of  dimensionality,  e.g.  in  canonical  variables 
(Anderson,  1958)  and  generalized  canonical  variables  analysis  (Sen  Gupta, 
1981b,c).  Finally,  if  the  means  and  variances  are  known  and  hence  can  be 
considered  as  zeros  and  ones  respectively,  the  discussions  in  the  previous 
paragraph  translates  to  the  case  of  SSMN  distributions.  The  present  case 
is  also  quite  interesting  from  several  theoretical  considerations.  Firstly, 
it  provides  a  practical  example  of  Efron's  curved  exponential  family 
(Efron  &  Hinkley,  1978)  and  illustrates  some  associated  difficulties  and 
techniques  to  overcome  them.  Secondly,  observe  that  this  situation  arises 
when  the  components  of  X  are  exchangeable,  i.e.  the  distribution  of  X 
is  independent  of  permutation  of  its  components.  Also,  the  SSMN  distribution 
constitutes  an  example  of  the  mean- zero  invariant  model  of  Andersson. 


Thirdly,  it  indicates  (through  p)  the  construction  of  simple,  intuitive 
estimators  of  the  correlation  coefficient  for  such  models  in  contrast 
to  those  obtained  by  the  iterative  method  of  scoring  of  Fisher  or  the 
sequential  Robbins-Monro  procedure  suggested  by  Sampson  (1976)-  Finally, 
it  is  demonstrated  that  unlike  the  LR  tests,  small  sample  optimal  test  for 
the  correlation  coefficient  may  be  conveniently  derived  in  such  models. 

Though,  as  noted  by  Anderson  (1963)  ’the  theory  in  the  case  of  correlation 
matrices  is  much  more  complicated  than  for  covariance  matrices  and  no 
general  result  could  be  given  in  a  simple  form',  we  observe  that  for  the 
above  important  and  special  structure  of  the  correlation  matrix,  interesting 
and  elegant  results  can  be  derived. 

Several  estimators  and  tests  for  p  are  considered.  Like  the  bivariate 
case,  the  maximum  liklihood  estimator  (MLE)  is  obtained  as  a  root  of  a  cubic 
equation  and  is  shown  to  lie  in  the  interval  restricted  by  p.  But,  unlike 
the  incomplete  bivariate  case  (Dahiya  &  Korwar,  1980)  it  may  be  necessary 
to  evaluate  all  real  roots  of  this  equation  which  makes  the  estimation  cumber¬ 
some  and  the  exact  distribution  of  the  estimator  nearly  intractable.  This 
also  renders  the  LR  test  statistic  computationally  quite  inconvenient. 

Further,  it  is  shown  that  the  LR  tests  for  one-sided  alternatives  are 
vacuous  with  positive  probability.  This  serious  drawback,  coupled  with 
the  lack  of  knowledge  regarding  its  small  sample  properties  and  computational 
inconvenience  motivates  the  development  of  alternative  tests.  A  test  based 
on  the  best  natural  unbiased  estimator  (BNUE)  of  p  is  proposed.  Also, 
a  LMP  test  for  p  is  derived  and  is  seen  to  co> ncide  with  that  based  on 
the  BNUE  of  p.  This  test  is  proved  to  be  unbiased.  The  exact  distribution 
of  this  test  statistic,  historically,  happens  to  be  a  problem  attempted  by 


Pearson  et  al  (1932, p. 341) .  The  exact  and  asymptotic  distributions  are 


derived  here.  Percentage  points  are  available.  The  case  of  constrained 
parameter  space  is  considered. 

The  generalized  variance,  the  determinant  of  the  dispersion  matrix, 
was  proposed  by  Wilks  (1932)  as  a  scalar  measure  of  multidimensional 
scatter.  It  has  enjoyed  applications  in  nearly  all  branches  of  applied 
statistics.  But,  tests  and  associated  exact  distributions  (SenGupta,  1981a) 
become  quite  complicated  in  the  general  case.  However,  for  the  SSMN  distri¬ 
bution  some  simplications  are  available.  The  LR  test,  which  again  is 
computationally  inconvenient,  is  derived.  Alternatively,  an  elegant  and 
simple  test  based  on  a  characterization  through  the  characteristic  roots 
of  the  correlation  matrix  is  proposed  by  introducing  the  concept  of  conditional 
characterisitic  roots. 

2.  Tests  for  equi-correlation  coefficient  p. 

Let  Y  ,  ...,Y  be  an  independent  sample  from  N  (0,5!  )  where,  letting  I  and  E 
—1  — m  k  —  p 

be  the  identity  matrix  and  the  matrix  with  all  elements  equal  to  unity  respectively, 

Ep  =  (1-p)  I+pE;  E"1  =  (1-p)  _1I-p  ( 1-p)  _1{  l+(k-l)p}_1E  «(c  ), 

Cii  =  {1+<k-2)p}/(l-P){l+<k-l)p},  c..  *  -p/(l-p) {1+ (k-1) p} ,  i  +  j. 

Hence  the  density  function  for  non-singular  Ep  can  be  written  as 

!  r  1  <r  y\)  <Iy.)2(-P) 

f(-;p)  *  (2iT)k/2"|E  |1/2  SXP  ~  ^  +  <1+<k-1)P)  ■**’■  (2'1) 

-0O  <  <  cof  i  =  -l/(k-l)  <  p  <  1. 

The  above  representation  is  particularly  useful  because  it  shows  that 
(i)  the  density  function  constitutes  a  member  of  Efron's  curved 
exponential  family 

(ii)  there  does  not  exist  any  one-dimensional  sufficient  statistic  for  p 


—  2  — 

(iii)  £  (yi~y)  and  y  *  £  y^/k  at®  independent,  and 
(iv)  the  part  of  the  exponent  within  the  second  bracket  is  monotonically 
decreasing  in  p,  with  positive  probability. 

We  want  to  test  HQ:  p  =  pQ  against  p<(>)pQ  or  against  H2:  p  f  pQ. 

2.1  Likelihood  ratio  test.  For  testing  against  the  LRT  is  derived 
below.  The  LRT  can  be  performed  even  if  ■  *  1.  The  likelihood  function  can 
be  written  from  (2.1)  above  easily  and  differentiating  this  with  respect 
to  p  and  equating  the  derivative  to  zero  we  have 

g (p)  *  (k-l)kp(l-p)  {l+(k-3)]p+££  y2 . {l+(k-l) p}2+£ (£  y.  ,)2{l+(k-l)p2}  =  0, 

ji  13  j  i  13 


where  =  (Y  ,...,y  ),  j  =  l,...,m.  This  is  a  cubic  equation  in  p, 

“~3  -13  *3 

two  of  whose  roots  may  be  complex.  Now,  g{-l/(k-l)},  g(l)  and  g(0)  are  all 
positive  (with  probability  one).  Thus,  it  is  not  clear  from  g(p)  that  the 
ML  estimator  of  p,  say  p,  with  '•l/ (k-lKjS'Sl ,  will  always  exist.  However,  it 
can  be  shown  directly  that  lim  f(y;p)=  0  for  p  {-l/(k-l)}+  or  p  +  1  which 
also  conforms  to  a  general  result  of  Anderson  (1970) .  Since  a  cubic  equation 
must  have  at  least  one  real  root  it  follows  then  that  the  ML  estimator  must 
be  in  the  interval  restricted  to  p.  Further,  in  cases  of  several  admissible 
solutions  (at  most  three) ,  by  principle  of  Maximum  Likelihood,  we  choose 
as  the  MLE  of  p,  that  which  corresponds  to  the  largest  value  of  the  Likelihood 

A 

function  and  call  it  p.  Thus  we  get  the  following. 

Theorem  1.  Let  Y  ~  N.  (0,£  ).  Then  if  Y, , . . . , Y  constitutes  an  independent 
- - -  —  k  —  p  — i  — m 

random  sample  from  the  above  population,  the  Likelihood  Ratio  Test  for 
testing  Hfl:  P“Pq  against  the  alternative  H2?  P^Pq»  is  given  by 

Reject  Hq  iff  X  «  {|£^|/|£p  | }m/,2exp[-  ~  {a{f ^ (P0)-f ^ (P) }+b{f2 (pQ) -f 2 (p) }}] 

<  K,  a  constant 


-6- 

where  p  is  the  MLE  of  p,  a  =  E  E  y2 , ,  b  =  E(E  y..)2,f^(p)  =  l/(l-p), 

j  i  13  j  i  1D 

f2(P)  =  -p/{l+(k-l) p} (1-p)  and  K  is  a  constant  to  be  determined  so  that 

the  level  of  the  test  meets  the  specified  value. 

2 

Under  H^,  for  large  m,  -2  lnX  ~  It  c^ear  that  the  ML  estimation 

A 

for  p  is  cumbersome  and  the  exact  distributions  of  p  and  the  LR  statistics 
are  nearly  intractable. 

For  one-sided  alternatives,  the  small-sample  behaviors  of  the  LR  tests 
can  be  even  worse.  For  example,  defining  a  test  to  be  vacuous  if  the  test 
statistic  is  a  constant,  we  have, 

Theorem  2.  The  likelihood  ratio  tests  for  testing  HQ^:  p  =  p+(0  <  p+  <  1) 
against  H^:  p  <  p+  and  also  for  testing  H^:  p  =  p  (-l/(k-l)  <  p  <  0) 
against  H12:  P  >  P  are  vacuous  with  positive  probability. 

Proof  Let  m  =  1,  since  similar  proof  holds  for  m  >  1.  Now, 

| £  |  *  (l+(k-l)p) (l-p)k_1,  d|Ep|/dp  =  -k(k-l) p (1-p) k-2 .  So,  | Zp |_1  t  for 
0  <  p  <  1  and  I  for  -l/(k-l)<  p  <  0.  Next,  consider  the  representation  of 
Y* as  *n  the  Potion  within  the  second  bracket  in  (2.1).  Then 

d(Y’E"1Y)/dp  =  [{l+(k-l) p}2 (Ey2)  +  {-1- (k-1) p2} (EyJ 2]/{ (l+(k-l) p) (1-p) }2 

=  [Cl(p)(Ey2)  +  c2(p) (Ey.)2]/{(l+(k-l)p) (1-p)}2 

2 

=  h(p,Eyi,Eyi) ,  say. 

Consider  testing  HQl  against  H^.  Observe  that  if  Y/EpSf  +  P  f°r  0  <  p  <  1, 
then  Sup  f(*)  under  will  be  attained  at  p  *  p+  which  also  yields  Sup  f(*) 
under  H^.  Hence,  to  prove  the  theorem  for  these  hypotheses,  it  suffices  to 
show  that,  h(»)  <  0  with  positive  probability . 


By  the  reduction  to  the  canonical  form  (Rao,  1973)  for  SSMN  distribution 


2  2 

there  exists  an  orthogonal  transformation,  Y_  -*■  Z,  such  that,  Zy.  =  Zz.  and 
z^  =  Zy^//k  ,  where  z^#  i=l,...,k  are  all  independent.  It  follows  that 
~  N(0,l+(k-l)p)  and  z^  ~  N(0,l-p),  j  =  2,...,k.  Hence 

p[h(p,Zy^,Zyi)  <  0]  =  Pftkc^pJ+c^p) }{l+(k-l) p>x^+(l-p)c1 (p)  <  0] 


2  2  2 

where  an<*  Xk  ^  are  independent  x  variables  with  1  and  k-1  degrees  of 

freedom  respectively.  To  show  that  the  above  probability  is  positive,  it 

2 

suffices  to  show  that,  kc^pj+c^p)  =  -(k-1)  (1-p)  is  negative  under  H  , 
which  is  obviously  true.  This  establishes  the  theorem  for  testing  Hq1  against 
H^.  Since  h(*)  is  an  indefinite  quadratic  form  in  Z^’s,  and  c^(p)  >  0  and 
c2(p)  <  0,  suitable  modificationsof  the  above  arguments  also  establish  the 
theorem  for  testing  HQ2  against  H^.  Due  to  the  above  difficulties  we  consider 
below  several  alterantive  tests  for  p. 

2.2  Test  based  on  BNUE  of  p.  From  (2.1),  where  m=l,  it  follows  that 
2 

(Zy^,ZyJ  is  a  sufficient  statistic  for  p.  But  this  is  not  complete,  since, 

E(Zy^)  =  0.  Note  that,  E(y^y^,)  =  p,  i  ^  i',  i,i*  =  l,...,k.  Further, 

2  2 

(Zy.)  -  Zy.  =  Zy.yJ,  so  that  based  on  a  sufficient  statistic  a'natural' 

l  l  l  i 

unbiased  estimator  for  p  is  Z  y^y^,/k(k-l) .  For  m  >  1,  considering 


i*i' 

n,  Z 
j 

the  best  (minimum  variance)  natural  unbiased  estimator  (BNUE)  of  p  is 


natural  estimators  of  the  form,  Z  a.(  Z  y. .y. . ./k(k-l) )  it  follows  that 

j  3  i*i'  13  13 


P  =  £  Z  y.  .y. ,  ./ink(k-l) . 
j  ijrfi'  3  3 

Hence  a  test  for  p  against  one  or  two-sided  alternatives  can  be  based 
on  p.  But  it  is  known  that  a  test  based  on  a  good  estimator  need  not  be 
a  good  test  [e.g.  a  test  for  correlation  coefficent  in  the  bivariate  case 
(Kendall  &  Stuart  (1967)].  However,  for  the  present  case,  it  is 

reassuring  to  note  the  following  desirable  result. 


-8- 


Theorem  3.  For  testing  H  :  p=0  against  H, :  p  <  (>)0,  the  test  which  rejects 
________  o  l 


H  iff,  p  <  (>  c')c,  where  c(c')  is  determined  to  give  the  desired  level 
o 


of  significance,  is  unbiased. 

Proof:  For  the  canonical  form  of  SSMN  distribution  discussed  in  Theorem  2, 


using  m  similar  orthogonal  transformations  on  y^'s,  i=l,...,k,  one  for 


each  j,  j=l,...,m,  Y.  -*■  Z  , ,  we  have, 
"3  “3 


p  -  I{fZ  y..)2  -  a  y2.)}/mk(k-l)  =  E{k  z2  -E.z2 .  }/mk(k-l) 
.  i  13  L  13  j  lj  i  13 


(2.2) 


~  2  2 
Then,  p  is  distributed  as  { (l+(k-l)p)  (k-Dy^  -  (l-p))^  }/mk(k-l) 


2  2  2 

where  X^  and  ^  are  independent  X  variables  with  m  and  m(k-l)  d.f. 


respectively.  Under  H  ,  the  distribution  of  p  is  the  same  as  above,  with 

o 


p=0. 


Consider  H^:  p  >  0.  To  prove  the  theorem  in  this  case,  it  suffices  to 


show  that 


which  is  clearly  true.  The  proof  for  :  p  <  0  follows  similarly. 

Note  that,  -  «  <  p  <  oo  and  it  may  be  desirable  to  consider  the  modified 
truncated  estimator  p,  where 


f-l/(k-l) 
P 

v.  1 


p  <  -l/(k-l) 
-l/(k-l)  <  p  <  1 
P  >  1 


-9- 


2.3  Locally  most  powerful  test.  No  small  sample  optimal  test  for  p 
seems  to  be  known.  One  such  test  is  presented  in  the  following  theorem. 
Theorem  4.  The  unbiased  tests,  defined  in  Theorem  3,  based  on  the  best 
natural  unbiased  estimator  p  of  p,  are  also  locally  most  powerful  tests. 
Proof :  Let  denote  the  first  derivative  of  the  log-likelihood  function 
with  respect  to  p.  Then  a  LMP  test  for  testing  HQ:  p=0  against 
H^:  p  >  (<)  0  is  given  by, 

SLq  >  (<)  k  ,  (2.3) 

where  k  is  to  be  determined  to  provide  the  desired  level  of  the  test. 

Now,  (2.3)  is  equivalent  to 

-{^(0)  +  ^(O)}  >  (<)  k' 

where  (p)  «  dll  |/dp  and  T-tp)  -  d(£  Y.E-1  Y.)/dp.  Use  of  the  expressions 
1  '  P  2  —j  p  —j 

for  ^(p)  and  ^(p)  as  given  in  proof  of  Theorem  2  and  some  simplifications 
establish  the  theorem. 

Due  to  the  remark  (i)  following  (2.1),  from  Kallenberg  (1981)  we  conclude 

that  the  shortcoming  of  the  LMP  test,  under  suitable  conditions,  tends  to  zero 

3/2 

at  the  rate  m  I  log  a  |  where  a  e(0,l)  is  the  level  of  significance. 

m  m 

2.4.  Exact  null  and  non-null  distribution  of  p.  The  exact  distribution  of 

2 

is  that  of  the  weighted  difference  of  two  independent  x  variables  with 
different  weights  and  arbitrary  d.f.s.  Now,  historically,  this  problem  was 
attempted  by  Pearson  et  al  (1932,  p.  341)  and  later  solved  only  partly  for 
the  very  special  case  of  equal  weights  and  equal  d.f.s  by  Pachares  (1952). 

It  was  also  encountered  by  Anderson  (1963,  p.  139)  who  conjectured  a  possible 
approximation.  The  distribution  is  presented  below  in  terms  of  Kummer’s 
function  .  For  percentage  points  see  Gpkhale  and  SenGupta  (1982) . 

Let  U(a,b; z)  give  independent  solutions  to  the  confluent  hypergeometric 
differential  equation  of  Kumnter: 

2  2 

zd  u/dz  +  (b-z)  du/dz  -  aa)  ■  0. 


■o< 


Then,  in  terms  of  an  integral 


U(a,b;z)  =  [<T (a)]  1  /  exp(-zv)ya  1(l+V)b_a_1  dv  ,  a  >  0,  z  >  0 

o 

and  in  terms  of  the  hypergeometric  function, 

IT  rlVa'b;z)  zl_blF1(1+a-b'2-b'2) 

'  sin  TTb  T(l+a-b)r(b)  "  r(a)T(2-b) 

Theorem  5 .  Let  V  =  Ct^V^  -  a2V2  w*iere  ai,012  >  ®  an<^  are  independent 

2  . 

X  variables  with  and  \>2  d.f.  respectively.  Then,  the  probability  density 
function  of  V  is  given  by, 

(V  +V  -2)/2 

f(v)  =  [c(v1,v2)/r(v1/2)]v  1  z  exp(-v/2a  ) 

•u[v2/2,  (V1+V2)/2;  {  (Ctj+O^ )  /2a  a2  }v]  ,  v>0 

(V  +v,-2)/2 

=  [c(vi,v2)/r(v2/2)](-v)  z  exp(v/2a2) 

•U^/2,  (v1+v2)/2;  {-<a1+a2)/2a1a2}v],  v  <  0 

(v  +V  )/2  V  /2  V  /2 

where  C  (v  ,v  )  =2  a.  az  . 

12  12 

00 

Proof ;  axf(v)  =  /  g1(  (v+a2v2)/a1>  g2(v2>  dv2 

where  g^  represents  the  probability  density  function  of  V  ,  i=l,2. 

For  V  0,  noting  that  the  limits  of  the  above  integral  reduce  to  0  and  », 
simplifications  yield  the  form  of  f(v)  as  in  the  theorem.  For  v  £  0, 
in  order  to  represent  f(v)  in  terms  of  U(*)  consider  the  following.  First 
note  that  the  exponent  for  U(*)  must  be  negative.  An  initial  transformation, 
<*2V2  =  -vy  and  then  a  further  tranformation  y-1  =  z  yields  the  claimed  result. 

Using  the  above  theorem  and  (2.2)  we  have  the  following  Corollary. 
Corollary .  The  exact  non-null  distribution  of  p  is  given  by  f(v)  of  Theorem  5 
with  =  m,  V2  =  m(k-l),  =  {l+(k-l)  p}/mk  and  a  =  (l-p)/mk(k-l)  from 

which  the  null  distribution  is  obtained  by  substituting  p  =  0. 
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2.5  Asymptotic  null  and  non-null  distributions  of  p  .  From  the  representa- 

~  2 
tion  of  p  in  Section  2.2  as  a  weighted  difference  of  two  independent  x 

variables,  it  follows  that 


E(ph)  =  l  (-l)h”3(^)  a3^-3^' 


j=0 


1  2  Hl, j  2,h-j ' 


h=l , 2, . . . 


•  2s  *  2s 

where  and  a2  are  given  in  Section  2.4,  y.^  E(X^  )  and  V2  g  *  E(x^k 

_  m  '  ' 

s  =  0,1,..., h.  Recalling  that  p  =  £  {  Z  y . .y . , ,/k (k-1) }/m,  by  Central 

j=l  i^i *  13  1  3 

limit  theorem,  we  have 

h  ~ 

Theorem  6.  k(m/2)  (p-p)  is  distributed  asymptotically  as  a  normal 

variable  with  mean  0  and  variance  {l+(k-])p}^  +  (1-p) 2/(k-l)  . 

2.6  Tests  for  p  under  constrained  parameter  space.  In  many  applied 
problems,  in  addition  to  having  SMN  or  SSMN  distribution,  further  informa¬ 
tion  on  p  may  be  available.  For  large  k,  we  need  p  >  0  in  order  that  Z^ 
be  positive  definite.  Again,  e.g.,  in  the  problem  of  psychological  testing 
theory  considered  by  Wilks  (1946)  or  in  testing  (Box,  1950)  the  Model  II 
assumptions  for  a  balanced  one-way  analysis  of  variance,  etc.  it  is  natural 
to  require  p  >  0.  Gleser  and  Olkin  (1969)  have  discussed  these  examples  and 
derived  likelihood  ratio  tests  for  p  under  the  restriction  p  >  0  on  the 
parameter  space.  We  consider  the  same  problem  for  SSMN  distributions.  As 
noted  earlier,  even  with  unconstrained  parameter  space,  the  likelihood  ratio 
test  suffers  from  several  undesirable  properties.  We  propose  below  an 
alternative  test  based  on  the  characteristic  roots  of  the  sample  correlation 
matrix,  study  its  unbiasedness  and  provide  the  null  and  non-null  asymptotic 
distributions  of  the  test  statistic.  Through  obvious  modifications,  the 
same  procedure  can  be  applied  to  the  case  of  p  <  0. 

Let  the  smallest  characteristic  root  of  Z  be  X.  Then  for  p  >  0, 

P 

X-l-p  with  multiplicity  k-1.  We  may  estimate  X  by  1-p  where  p  is  a 


-12- 


modified  estimator  of  p,  e.g.  the  modified  maximum  likelihood  estimator  or  the 

modified  BNUE,  modified  by  the  restriction  p  >  0.  Alternatively  consider 

k 

estimating  X  directly.  Now  d  =  £  d./(k-l)  is  an  intuitive  estimator 

i=2  1 

of  X,  where  d  >.  • • •  >  d  are  the  characteristic  roots  of  the  sample 
correlation  matrix.  Further,  from  Anderson  (1963),  [ (m-1) ** (d-X)/{X (k- (k-1) X}] 
[k(k-l)/2]^  is  asymptotically  a  standard  normal  variate. 

Theorem  7.  Let  Y_,...,Y  be  a  random  sample  from  the  SSMN  distribution 
-  — 1  — m 

with  p  >  0.  Then 

(a)  Testing  HQ:  p=pQ  against  H  :  p  >  (<)  pQ  is  equivalent  to  testing 

H0:  ^  =  Xo=1-po  ai3ainst  H{:  ^  <  <>)  ^0- 

(b)  The  test  which  rejects  HQ  in  favor  of  p  >  (<)p  if  and  only 

if  d  <(>  d^)  dp f d p ( d p )  a  suitable  constant,  is  asymptotically  un¬ 
biased  for  all  p  satisfying  0  <  pQ  <  p  <  (k-2)/2(k-l) 

[(k-2)/2(k-l)  <  p  <  p0  <  1]. 

Proof :(a)  First  note  that  for  a  k  x  k  correlation  matrix,  (k-1)  roots  are 
all  equal  if  and  only  if  it  is  of  the  form  The  if  part  is 

trivial.  The  only  if  part  follows  from  Anderson  (1963,  Appendix 
A).  The  desired  equivalence  then  follows  trivially. 

(b)  Consider  H^:  p  >  PQ  >  0.  We  use  the  asymptotic  distribution  of  d 
given  above.  Then,  for  large  m, 

P(d  <  d0|Hp)  =  a  <=>  dp  =  (-T  )  Xp  (k-  (k-1)  XQ)  [2/k  (k-1)  (m-1)  ]**  +  XQ 
The  proposed  test  will  be  asymptotically  unbiased  if,  for  large  m, 

P (d  <  dp | H ^ )  >  a,  i.e.  if,  for  all  X^  <  XQ, 

(-T  )  Xp{k-  (k-D  Xp}[2/k(k-l)  (m-D  l*5  +  (Aq-X^ 

>  C-T  XX  fk- (k-1) X  } .  [2/k (k-1)  (m-1)]15 

a  1  1 

which  reduces  to  the  condition  Xp{k- (k-1)  XQ}  <  X^k-  (k-1)  X^}  .  Let 

g(X)  »  X{k- (k-1) X} .  Then  g(X)  I  X  if  X  >  k/2(k-l)  or  p  <  (k-2) /2 (k-1) ,  as  claimed. 

The  proof  for  H| :  0  <  p  <  pQ  follows  as  in  above  with  obvious  modifications. 
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3.  Tests  for  the  generalized  variance  |£p|>  As  stated  in  Section  1  and 
also  discussed  by  Eaton  (1967)  tests  for  generalized  variance  are  often 
of  practical  importance  and  interest.  We  consider  several  tests  for  | | . 

3.1  Likelihood  ratio  test.  First  note  that,  | |  <  1  for  p  ?  0.  Next, 
consider  the  following 

Lemma  1.  If  |l  |  =  cf^  <  1,  then  there  are  precisely  two  real  distinct 
solutions  for  P. 

Proof;  Recall  (from  the  proof  of  Theorem  2)  that  |E^|  t  for  -l/(k-l)  <  p  <  0 
and  +  for  0  <  p  <  1.  Hence,  due  to  strict  monotonicity  there  are  two  real 
solutions,  say  p2  >  0  >  p^  to  |E  |  =  aQ. 


Then  (using  MLE  of  p  from  Section  2.1)  we  have 


2  2 

Theorem  8 .  The  likelihood  ratio  test  for  H_:  |E  I  ■  a  against  h  :  |E  |  ?  o'! 
-  0  1  p  0  1  p  C 


is  given  by 


Reject  H0  iff  Xm  =  ( |  l/ao)m/2  exP{_Js  l  c 


where  p**  is  such  that  II  f(Y.?p**)  =»  max  IT  f  (Y . '•  p) ,  p  is  the  MLE  of  p  and 

1  P1'P2  1 

C  is  a  constant  to  be  determined  such  that  the  test  has  the  desired  level. 

2 

Under  HQ,  for  large  m,  -2  In  Xm  'VX^* 

The  LR  test  statistic  is  computationally  cumbersome.  For  one-sided 
alternatives,  it  becomes  even  more  repulsive.  Alternatively,  another  test 
is  considered  below. 

3.2  Tests  based  on  the  smallest  characteristic  root.  Let  p^  and  P2  be  as 

defined  in  Proof  of  Lemma  1.  Now,  the  distinct  characteristic  roots  (c.r.) 

for  E  are  l-pj  and  l+(k-l)p. ,  i  -  1,2.  Thus  the  following  definition  is 
pi  1  1 

introduced,  noting  that  for  (k-l)pi  ?  -p  ,  the  c.r.'s  are  all  different. 


-14- 


Definition.  1-p^,  l+(k-l)p^,  i=l,2  are  defined  as  the  conditional 

2 

characteristic  roots  (c.c.r)  of  Z  subject  to  I Z  I  =  a  ,  where  p.,  i=l,2 

p  1  p1  0  l 

2 

are  the  two  distinct  real  solutions  to  I Z  I  =  CL  . 

1  p  0 

2 

Lemma  2 .  The  c.c.r.s.  are  all  different  for  0<a  <1,  k>2. 

Proof ;  It  suffices  to  check  C^:  (k-Dp^  ?  -p ^ ^  j'i*}"!*2* 

Consider  C^.  Suppose,  if  possible,  (k-ljp^  =  “Pj-  Then,  from  Lemma  1, 
the  defining  relation  between  p^  and  p^,  i.e.  jZ^  |  =  |Z^  j  reduces  to 

T(P2)  =  [{l+P2/(k-l)}/U-p2)]k-1  -  [{l+(k-l)p2}/(l-p2)]  =  0. 

To  establish  it  suffices  to  show  that  ¥(p)  has  no  real  solution  in 
(0,1).  Now,  letting 


ftp)  =  f*  <p)  -  f2(P>.  fx(P)  =  {l+p/(k-l)}/(l-p),  f2(p)  =  {l+(k-l)p}/(l-p) 


we  have 

V (P)  =  [fk_2(p)-l]k/(l-p)2. 

But,  p  >  0  <=«>  f^p)  >  1  <=>  (fk_2(p)  >  1,  k  >  2}  <=>  f’(p)  >  0, 
so  that  f  is  strictly  monotonically  increasing  for  0  <  p  <  1.  The  desired 
result  then  follows  by  noting  that  f(p)  is  continuous  at  p  *  0  with  f(0)  *  0. 
A  similar  argument  yields  C  since  0  <  f  (p)  <1  for  -l/(k-l)  <  p  <  0. 

2 

For  k=2,  p^  =  -p2  so  that,  testing  ) Z^ |  is  equivalent  to  testing  p  . 

Thus,  in  the  following  sequel,  it  is  assumed  that  k  >  2.  Then,  the  smallest 

2  2 

c.c.r.,  s.c.c.r(Zp| |Zp|  *  aQ  ,aQ  known)  is  uniquely  defined  by  virtue  of 
Lemma  2  and  it  will  be  denoted  by  X(p,aQ). 


Lemtna  3»  If  lEyl  *  ax2  <  aQ2  =  |Ep|  <  1,  then  Xjy,^)  <  X(p,0Q). 

2 

Proof:  First  consider,  |Z  |  *  a_  and  let  o.  <  P-  be  the  two  real 

-  ■  ■  ■  ■  p  0  12 

solutions.  l+Oc-Dpj^  and  l-p2  are  less  than  both  1-p^^  and  l4(k-l)P2  since 
P1  <  0  <  p2>  Then  the  s.c.c.r  is  either  l+Oc-lJp^  or  l-p2> 

Next,  observe  that  X(p,ap)  and  X(y ta^)  must  be  of  the  same  functional 
form.  Consider,  e.g.,  the  case,  X(p,Cp)  **  l+Ot-Dp^,  i.e.  (k-Dp^c-p^ 

If  possible,  suppose  Xfy,^)  -  1-V2'  i*e*  >  -Y2-  Here  k  is  fixed. 

Since  X(S,a  ),  -l/(k-l)  <  6  <  1,  is  a  continuous  function  of  a  ,  for 
(k-l)y1  >  -y2  to  hold,  there  must  exist  a  n  such  that  (k-l)n1“-Ti2, 
which  contradicts  Lemma  2. 

It  remains  to  show  that  l+tk-l^  >  1  +(k-l)y1  and  l-p2  >  l-Yj. 

This  follows  from  the  strict  monotonicity  of  |E^|. 

Theorem  9.  Testing  Hq:  | |  «*  oQ  against  :  |  |  «=  >  oq  is 

equivalent  to  testing  H^:  X(p,<Tp)  <  X(p,0^)  . 

Proof ;  For  a  k  x  k  correlation  matix,  (k-1)  roots  are  all  equal  iff 
it  is  of  the  form  Ep.  The  if  part  is  trivial.  The  only  if  part  follows 
from  Anderson  (1963,  Appendix  A)  as  also  stated  by  Lawley  (1963).  The 
theorem  then  follows  by  noting  Lemma  3. 

The  above  theorem  motivates  tests  for  |Z  I  to  be  based  on  estimators 
of  smallest  characteristic  root  (s.c.r.).  We  propose  the  tests: 

Reject  Hq:  | Ep|  -  0Q2  against  1^:  |Ep|  <  (>)  0Q2  iff  X  <(>)X0, 

and  against  H£:  |Ep|  i*  cp2  iff  X  <  X'  or  >  X" 


where  X  is  a  suitable  estimator  of  s.c.r.  and  XQ,X'  and  X"  are  chosen  so 
as  to  give  the  desired  level. 

In  the  following  sequel,  Lemmas  1,2  and  3  are  exploited  to  provide 
some  simple  tests  for  |Ep|.  Let  ^  *  M*-n* (1+ (k-1) P» 1-p)  where  p 

*y  A  _ 

is  a  suitable  estimator  of  4  e,g.  p,por-p.  Consider  p  =  p  the  truncated  BNUE 
of  p,  due  to  the  difficulties  associated  with  MLE  as  seen  in  Section  2. 

Then,  for  0  <  X  <  1, 

P  (X  <  XQ)  =  p{p  <  (X0-l)/(k-l)  |p  <  0}  P(p  <  0)  +  p{p  >  l-X  |p  >  0}  P(p  >  0). 

The  null  and  non-null  exact  distributions  can  then  be  obtained  from  the 
distribution  of  p  given  in  Section  2.4,  since  from  Lemma  2  the  s.c.c.r.  and  hence  p 
is  uniquely  specified  under  HQ  and  H^.  The  unbiasedness  of  this  test  can 
be  investigated  in  the  same  manner  as  in  Section  2.3. 

Hq  may  motivate  modified  estimators  of  s.c.r.  Suppose,  e.g.,  HQ 
specifies  X(p,oQ)  =  1  +(k-l)p^.  In  this  case,  let  us  assume  the  population 
s.c.r.  is  of  the  form  l+(k-l)p,  p  <  0.  Since,  also  X(p,0Q)  and  XCYjO^) 
must  be  of  the  same  form  as  proved  in  Lemma  3,  one  may  restrict  to  estimators 

of  the  form  l+(k-l)p.  One  approach  could  be  to  use  the  modified  BNUE  or  MLE  of  p, 
modified  by  the  restriction  p  <  0.  Another  would  be  to  consider  the 

characteristic  roots  of  the  sample  correlation  matrix  R.  For  example, 
being  the  smallest  root  also  specifies  that  the  s.c.r.  in  the 
population  is  of  multiplicity  1.  This  enables  us  to  use  Anderson's  (1963) 
result:  [  (m-1)  ** (d-Xj) /{X2  (k-q^)  >]  <k.q1<q2/2)  **  is  asymptotically  distributed 

as  a  standard  normal  variable  where  the  k  x  k  population  correlation  matrix 
R  has  q^  roots  equal  to  Xj  and  q2  *  k-q1  roots  equal  to  X2  <  X^,  d  ■ 

^  di'/,q2'  dl  —  '**  —  dk  are  characteristic  roots  of  R  based  on  a  sample 
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o£  size  m.  Hence,  for  large  samples,  the  test  can  be  based  on  d.  For 

A 

X(p,aQ)  ■  l+(k-l)p1  or  1“P2#  q2  -  1  or  k-1  respectively.  Further,  A 

can  be  taken  as  d  and  a  test  for  H„s  IE  I  *  0.  can  be  based  on 

0  P  0 

it  as  before.  The  asymptotic  distribution  of  d  is  completely  specified 
under  HQ  and  whenever  the  corresponding  s.c.c.r.s  are  defined. 

4.  Applications .  The  model  in  this  paper  assumed  known  means  and  variances. 
Wilks  (1946)  considered  the  same  model  with  known  means  and  unknown  variances. 
Cox  (1958)  gives  a  practical  example  of  known  mean  and  known  variance  situa¬ 
tion.  Efron  and  Hinkley  (1978)  discuss  the  same  model  as  in  this  paper,  for 
the  special  case  of  k  =  2  in  some  detail  in  terms  of  information  and  curvature. 
Goodman  (1981)  applies  this  model,  with  k»2,  to  the  analysis  of  two-way 
contingency  tables.  The  SSMN  distribution  is  also  often  applied  to  problems 
in  reliability.  The  methods  discussed  here  are  hoped  to  also  indicate  solutions 
to  a  variety  of  problems  with  patterned  correlation  matrices,  e.g.  multi¬ 
variate  intraclass  correlation,  autoregressive,  moving  average  etc.  models 
(Sampson  (1978)).  For  detailed  discussions  and  a  new  development  using  the 
notion  of  mean  curvature,  of  the  multiparameter  LMP  test  see  Sen  Gupta  and 
Vermeire  (1981).  This  may  be  used  to  generalize  the  results  of  Section  2.3 
to  the  multivariate  intraclass  correlation  models. 
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